CP violation for neutral charmed meson decays to CP eigenstates^] 
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CP asymmetries for neutral charmed meson decays to CP eigenstates are carefully studied. The 
formulas and numerical results are presented. The impact on experiments is briefly discussed. 
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Up to now, we still do not have any experimental evidence for CP violation in the charm sector. Theoretically, 
the prediction for charm mixing is very small. This leads to small CP violating effects in charm decays. However, 
searching for large mixing and CP violation in charm decays is still very interesting not only for testing the standard 
model but also for finding new physics( for a recent review, see Ref.[l|). Because CP eigenstates are very special, if 
D° — D° decay into the same CP eigenstates, then the CP violating asymmetry could be enhanced by interference. 
So we try to investigate this possibility in detail. Another advantage of CP eigenstates is that the amplitude ratio 
A(D° — > f)/A(D° — > /) can be estimated without computing the amplitudes directly. This makes the computation 
of the CP asymmetries easier. In this paper, we shall concentrate on the case of CP eigenstates into which charm 
decays. 



Define 



II. TIME-DEPENDENT CP ASYMMETRY 



CP \D°) = \W) , 
\D s )=p\D°)+q\W} , 

\D L )= P \D°)- q \W}, 

\p\ 2 + \q\ 2 = l. 



(1) 



The corresponding eigenvalues of \D S ), \D L ) are 



Ac = TOo — I— . At = TOt — I 



X 



Assuming CPT invariance, the time-evolved states are 



\D p (t))= g+ (t)\D ) + lg_(t)\D°) 



where 



\D°(t)) = Zg_(t)\D°)+g + (t)\D°) 



(2) 



— _ p—vmt—^t I pi^P-t- 

~ 2 



Am = to l — TOg , to = (m L + TO s )/2 , 
, A 7 = 7 S - 7l , 7 = ( 7l + 7s )/2 . 



(3) 
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Define the mixing parameter 



Am A7 
X = — < y= 2^< W 



then the decay amplitudes for final state / are 

A(D° p (t)^f) = (f\H eS \D° p (t)) 

= g + (t)A(f) + lg_(t)A(f) 

= A(f){g + (t) + X f g_(t)} , (5) 

where 

A(f) = (f\H cS \D°) , 
A(f) = (f\H eS \D°) , 

1 p A(f) 

Similarly, we put 

A(J) = (7_\H eS \D°) , 

A(f ) = (7\H e e\D°) , 

X-- P ^H (7) 

Xf ~ * 3(7) • (7) 

where / is the CP conjugate state of the final state /, and 

17) = CP I/) =r?cp (/)|/} , 

with ?7cp(/) = ±1 is the CP eigenvalue (or CP parity). For the amplitude of D^(t) — > /, we have from Eq. (2) 

A(W p (t)^J) = (J\H eS \W p (t)) 

= Zg-(t)A(f)+g + (tjA(f) 

= A(J){g + (t) + Xjg_(t)} . (8) 

Now, it is easy to calculate the time-dependent width T(D p (t) — > /) and T(D^(t) — > /). Using Eqs. (3), (5) and (8), 
we have 

T(D° p (t)^f) = \A(Dl(t)^f)\ 2 

= \A(f)\ 2 {\g + (t)\ 2 + 2Be[X f g* + (t)g_(t)] + \X f f\g_(t)\ 2 } , (9) 

r(D*(t)-7) = |A(7)| 2 {| ff+ (t)| 2 + 2Rc[Aj 5 ;(%_(t)] + |A 7 | 2 | ff _(t)| 2 } . (10) 

In order to compute Xf and Aj, we need first to compute the amplitude ratios A(f)/A(f) and A(f)/A(f). As 
an example, we consider D°,D° — > K + K~ . Draw the decay diagrams (Fig. 1), we see that if we neglect the 
penguin diagram contribution, the D° and D° decay diagrams involve only one CKM factor V US V* S and V* S V CS , 
respectively. The only difference of D° and D° decay diagrams is that the initial and final particles change into their 
CP counterparts. So 

A(f) A(W> -> K+K-) 
A(f) ~ A(D a — > K+K~) 

= VcP (K + K-)^^ = Vcp(K + K-) = +1 . (11) 

' 11 a 'ft 
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FIG. 1: Decay diagrams for D° , D° -> K + K~. 

In Eq. (11), V cs and V us are both real in Wolfenstein parametrization for CKM matrix and rj CP (f) is the CP parity 
of the final state /. Usually i]cp{f) — ±1 for different /. Actually, we can prove that (see the appendix in Ref.[2|) if 
the decays of D° and D° only involve one CKM factor respectively, then the ratio 



A(f) _ m £-^-_ m 
- VcpU) ~ VcpU) 



(12) 



The last equality holds only for charm decay because all the CKM matrix elements involved are real, if we neglect 
the penguin contribution. 



Define 



From Eqs. (6) and (7), we have 



Pf 



Mf) 

A(f) 



pj- 



A{f) 



(13) 



A/ - -P f =VcAf) 
Aj = ~Pj = Vcp(f) 



-lip 



(14) 



After a straightforward calculation we arrive at 



r(Dj!(i)-/) = i e -^|A(/)| 2 i(i + 



-|A 7 t j. e iA 7 t 



( ) + 2(l 



) cos Amt 



4 



+2r? C p(/) 



l e -^\A(f)\ 2 Ul + 



iA 7 t _ „|A 7 t 



(e-t^Tt _ e i^7t) C os 93 + 2 sin p sin Ami 
2 



(15) 



)( e -5^ t + e 2^) + 2 (l 



I cos Ami 



-2»7cp(/) 



e 2 7 ) cos ip — 2 sin ^ sin Amt 



Assume 



\A(f)\ = \A(f)\ 



(16) 



(17) 



This is guaranteed by our approximation of neglecting the penguin, because in that case only one CKM factor appears 
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The time-dependent CP asymmetry is 

r(D°(t)^f)-r(W p (t)^J)_N f (t) 



C f (t) 
N f (t) 



r(D°(t)-/) + r(D°(t)-/) D f (t) ' 



(18) 
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i cos Amt 
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■|A 7 t _ iA 7 t 



)cos^ + 4?7 C p(/)( 
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) sin ip sin Amt 



(19) 



) sin if sin Amt 



(20) 



In Eqs. (19) and (20), there are several parameters we need to know: the phase ip, x = Am/7, V — A7/27, M/|p|, 
etc. But we do know that \x\ ~ |y| < 10 _2 (Ref.[l[), and \q\/\p\ is very close to unity. Some people assume [J that 
M/IpI — H/M ^ ±1%- As for the phase ip, 



(21) 
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In B° — B° system, the box diagram dominance leads to 



q 

PJb 



-2i$ 



, pb = 2(3 



(22) 



In charm case, if we can use the same approximation and assume 6-quark plays the dominant role in the corresponding 
box diagram, then 



P/d 



-2i 7 



PD = 27 . 



(23) 



But it is not the case for charm. Firstly, the box diagram does not dominate. Secondly, even in the box diagram 
because \V ub \ is very small compared with \V ud \ and |V^ S |, the internal 6-quark contribution may be not important. 
Furthermore, V ud and V us do not carry the weak phase, cp can also get the contribution from the iT 12 /2 term. 
Anyway, we do not know the value of ip, so just keep it as a free parameter. For e ±2 A7 *, using Aj/(2j) — y, we have 
e ±iA 7 t = e ± vl t = e ±a*/r DOj because y < 10~ 2 , e ±yt/r v a is around unity. 



N f (t) c -4r Jcp (/)( 
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) (yjt) cos (p + 8?y C p (/) sin (p sin Amt , 



(24) 
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D f (t) 
C f (t) 



Vcp(f) 



Vcp(f) i 2^*( 
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P 



cos + sin tp sin Ami 
) cos tp + sin y> sin(a;7t) 



(25) 



(26) 



III. TIME-INTEGRATED CP ASYMMETRY 



In order to have more statistics, we integrate the time-dependent observables with time. We first list some useful 
quantities: 



G + = 



G 



G+- = 



a+ i u\V 2 + x 2 -y 2 1 

„ d * lfl+(t)l = 27(1 + ^(1-^)" 7' 



x 2 + y 2 



x 2 + y 2 



„ d * lff - ( * )|2 ~2 7 (l + ^)(l-y2)~ 27 ' 

ji * /.s /.\ -2/(1 + x 2 ) + zx(l-j/ 2 ) -j/ + mc 

dt g + (t)g_(t) = — — — — — - 

2 7 (l + x 2 ){l -y 2 ) 



2" 



for x , y <C 1. It is straightforward to get the integrated decay width. From Eqs. (9) and (10) we have 
T(D° p - /) = / di r(D°(t) - /) = |A(/)| 2 {G + + 2Re(A / G + _) + |A / | 2 G_} , 

JO 

T(Dl - /) = / df r(^(t) -> J) = L4(7)| 2 {G + + 2Re(A 7 G + _) + |Ay| 2 G_} . 
Jo 

Again we assume |A(/)| = |A(/)|, then the time-integrated CP asymmetry is 



r(Dg->/)-r(^-./) _ N f 

T(D^f)+T(D0^J) D } 



y( 



) cos tp — x( 



N f = -2?7cp(/) 

D f = 4-2ry CP (/) 
w 4 . 

Finally neglecting the (x 2 + y 2 ) term in Eq. (33), one can obtain 



) sin tp 
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P 




q 



) cos tp + x sin y> 



(27) 
(28) 
(29) 

(30) 
(31) 

(32) 
(33) 

(34) 

(35) 



In Ref. 3], the first term in Eq. (35) is omitted and the CP paritly factor ?ycp(/) is missing. 

Up to now, we have only discussed incoherent D° — D° decays. Sometimes D° — D° pairs are produced coherently 
such as in e + e~ colliding machines (BES and CLEO-c). 

The time-evolved coherent state of D° — D° pair can be written as [3] 



\i) = \D°(k 1 ,t 1 )D°(k 2 ,t 2 )) +r,\D°(k 1 ,t 1 )D°(k 2l t 2 )) , 



(36) 



where rj is the charge conjugation parity or orbital angular momentum parity of the D° — D° pair. 

Because D° — > l + X and D° — > l~X only, we can use the semileptonic decay to tagg one of the two time-evolved 
states D®(t) and D®(t). We define the leptonic-tagging CP asymmetry C« as 



C 



fi 



N(l-,f)-N(l+,f) 
N(l-,f)+N(l+J) 



(37) 
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where 



N(r,f) 



dt.dt, \(rj\H cS \z)f 



(38) 



is proportional to the number of events in which D®(k, t) — > I A as tagging in one side and the other side is the decay 
£>p(fc, t) — ► f or vice versa. Similarly, 



N(l+,f)= dt^Kl+JlH^ 



Assuming = |A(/)| and |A(Z + )| = \ A(l )|, after a tedious calculation, we have 

\A(l-)A{f)\ 2 {G\ + G 2 _ + 2\X f \ 2 [G + G_ + r,\G + _\ 2 ] 
+2(1 + r])G_Re(X f G* + _) + 2(1 + v )G + Re(\ f G + _) 
\A(l+)A(f)\ 2 {G\ + G 2 _ + 2\Jj\ 2 [G + G_ + rj\G + _ \ 2 
+2(1 + r 1 )G_Re(XjG* + _) + 2(1 + r 1 )G + Rc(XjG + _) 



N(r,f) 
N(l + J) 



2 V Re(G 2 + _)} 
2 V Re(G 2 + _)} 



Cfi - 
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fl 
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fl 



N(l-,f)-N(l+J) _ Nfi 
N(l-,f)+N(l+,f) D fl 
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(39) 



(40) 

(41) 
(42) 



(43) 



(44) 



Finally 



C 



N 



fl 



fl 



N 



(1 + v)VcpV) 



fl 



) cos ip + x sin ip 



(45) 



Comparing Eq. (45) with Eq. (35), we find that C« is just twice as large as C ^ when the charge conjugation parity 
or the orbital angular momentum I is even. This is surprising. From Eq. (35), the order of magnitude of Ci is 
< 1CP 3 , because x ~ y < 10~ 2 . Now we present C ^ (theory), C ^ (exp.), branching fractions for D°,D° decay into CP 
eigenstates and the number of D — D pairs needed for testing CP asymmetry for ltr signal lower bound in TABLE 
I, where for the branching ratios we take most of them from the 2006 particle data booklet [f|. For K* + K*~ , p + p~ 
and p°p° final states, we use the BSW theoretical estimation [6j. For the measured CP asymmetries listed in Table I 
are also taken from the 2006 particle data booklet We use the formula for N 



DD 



BC 2 



BC 2 



for lrr signature ; 
for 3cr signature . 
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TABLE I: The number of DD pairs needed for testing CP asymmetry 



D° ->fC f (theory) C f (exp.) BR N Q - (la lower bound) 

K + K~ 0.014 ±0.010 (3.84 ± 0.10) x 10~ 3 

K S K S -0.23 ±0.19 (3.7 ±0.7) x 10" 4 2.70 x 10 8 

K*+K*~ 1.0 x 10~ 2 (BSW) 1 x 10 7 

tt+tt" < 10" 3 0.013 ±0.012 (1.364 ± 0.032) x 10 -3 7.4 x 10 7 

7r°7r° 0.00 ±0.05 (7.9 ±0.8) x 10" 4 1.26 x 10 8 

p°7T (3.2 ± 0.4) x 10~ 3 3.13 x 10 7 

p + p- 1.3 x 10~ 2 (BSW) 7.69 x 10 6 

p°p° 1.2 x 10~ 3 (BSW) 8.33 x 10 7 

(j>ir° (7.4 ±0.5) x 10~ 4 1.35 x10 s 

4>rj (1.4 ± 0.4) x 10~ 4 7.14 x 10 8 

K*°K*° (7 ±5) x 10~ 5 



IV. SUMMARY AND CONCLUSION 



We have computed the time-dependent and time-integrated CP asymmetry for neutral charmed meson decays into 
CP eigenstates. We present CP asymmetry not only for incoherent D° — D°, but also coherent D°D° pairs. We find 
that the time-integrated CP asymmetries are very small (order of < 10~ 3 ). We also give the lower bound for the 
number of DD pairs needed for testing the CP asymmetries. At present, the integrated luminosities for e + e~ colliders 
are : 



BES II : 27 pb~ J 
BES III : 20 ft" 1 
CLEO - c : 281 pb~ 



for 4 years data taking 



The corresponding DD pairs are 



BES II : 
BES III : 
CLEO - c 



10 5 
10 7 
10 6 



In Table I, for D — ► VV decays, only when both vector mesons are longitudinally polarized the VV final states are 
CP eigen states. For the corresponding branching ratios in Table I, we assume that the VV final states for which 
both V are longitudinally polarized doninate. From Table I we see that the only hope is relying on BES III and 
B factories. At i?-factories, because the large data sample of charmed meson, both time-dependent asymmetry and 
time-integrated asymmetry can be measured. While for BES III, only time-integrated CP asymmetry can be tested. 
Of course, if there is new physics, some surprise may happen. We can also see from Table I that all the measured CP 
asymmetries are consistent with zero. 
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